Ab-initio approach to effective single-particle energies in doubly closed shell nuclei 
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The present work discusses, from an ab-initio standpoint, the definition, the meaning, and the 
usefulness of effective single-particle energies (ESPEs) in doubly closed shell nuclei. We perform 
coupled-cluster calculations to quantify to what extent selected closed-shell nuclei in the oxygen and 
calcium isotopic chains can effectively be mapped onto an effective independent-particle picture. 
To do so, we revisit in detail the notion of ESPEs in the context of strongly correlated many- 
nucleon systems and illustrate the necessity to extract ESPEs through the diagonalization of the 
centroid matrix, as originally argued by Baranger. For the purpose of illustration, we analyse the 
impact of correlations on observable one-nucleon separation energies and non-observable ESPEs in 
selected closed-shell oxygen and calcium isotopes. We then state and illustrate the non-observability 
of ESPEs. Similarly to spectroscopic factors, ESPEs can indeed be modified by a redefinition 
of inaccessible quantities while leaving actual observables unchanged. This leads to the absolute 
necessity to employ consistent structure and reaction models based on the same nuclear Hamiltonian 
to extract the shell structure in a meaningful fashion from experimental data. 



I. INTRODUCTION 

The concept of single-nucleon shells dates back to the 
early days of nuclear physics and constitutes the basic 
pillar of the nuclear shell model The independent- 
particle approximation provides a zeroth-order picture of 
the structure of nuclei on top of which correlations arc 
added to provide a more realistic description. Based on 
such a rationale, the correlated shell model has been able 
to explain the occurrence of extraordinarily stable config- 
urations for specific neutron and proton numbers, known 
as magic numbers. As a matter of fact, the universal 
character of such magic numbers over the nuclear chart 
remains an open question today . Recently, the evolu- 
tion of shell structure and the understanding of the neu- 
tron drip-line location in oxygen isotopes have received 
considerable experimental and theoretical attention 0- 
m, e.g. significant shell gaps have been identified in ^^O 
and ^^O leading to the interpretation of new magic shell 
closures at iV = 14, 16 in Z = 8 nuclei. 

Identifying the underlying mechanisms responsible for 
the occurrence or the disappearance of magic numbers 
in specific regions of the nuclear chart requires improve- 
ment on the traditional shell model by allowing for a 
more systematic and consistent inclusion of correlations. 
In particular, questions related to the impact of contin- 
uum degrees of freedom and of thrcc-nuclcon forces 
on the evolution of nuclear shells is a frontier driving 
low-energy nuclear physics research in connection with 
radioactive ion beam facilities [l3| • 



Whether a certain nucleon number qualifies as a (new) 
magic number cannot be postulated a priori. Experi- 
mentally, several quantities, e.g. the excitation energy 
and the collective character of the first 2+ state in even- 
even isotopes, the size of the gap in the one-nucleon addi- 
tion/removal spectrum, and the spectroscopic factors of 
associated low-lying states in odd-even neighbours need 
to be extracted in order to make such an assessment. 
Theoretically, the same quantities need to be computed 
while including all many-body correlations that could 
play a role in order to check whether the picture asso- 
ciated with a magic number eventually holds. 

It can be useful in this context, for analysis and in- 
terpretation purposes, to extract a single-nucleon shell 
structure, i.e. a set of effective single-particle ener- 
gies (ESPEs) associated with an underlying independent- 
particle-like picture the system is mapped on. How- 
ever, immediate non-trivial questions arise that arc at 
the heart of the present study 

1. Can a single-nucleon shell structure be unambigu- 
ously defined in a system that is intrinsically corre- 
lated? In other words, can ESPEs be computed 
on the sole basis of outputs of the many-body 
Schroedinger equation and not as a result of an a 
priori given zeroth-order approximation picture? 

2. Correspondingly, to which auxiliary independent- 
particle problem arc ESPEs related, i.e. which one- 
body Hamiltonian arc ESPEs the eigenvalues of? 

3. To what extent do correlations impact the effec- 
tive independent-particle picture provided by such 
ESPEs? 
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4. In which way are ESPEs related to underlying nu- 
clear forces? 
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5. Given that an unambiguous definition of ESPEs ex- 
ists, is tlie associated simplified picture needed and 
beneficial or potentially misleading? In particular, 
are ESPEs physically observable quantities? 

Several of the above questions have been answered 
long ago while others still necessitate further clarifica- 
tions. The procedure to extract ESPEs unambiguously 
(cf . point 1 above) goes back to French and Baranger [lll - 
|l3[ and can be utilized to address points 2, 3, 4, and 
5. Such a procedure defines ESPEs as centroid energies 
denoting barycenters of correlated total binding energy 
differences between the A-nucleon state the one-nucleon 
transfer is performed on and the complete set of eigen- 
states of the A-l-1 and A-1 systems. Eventually, cen- 
troid energies can be related [l^ - fl^ to the monopole 
part [m. Ha] of underlying nuclear interactions, which ef- 
fectively answers points 2 and 4 above. 

In spite of the existence of an unambiguous procedure 
to compute ESPEs, difficulties exist that can lead to im- 
proper conclusions, e.g. conclusions based on an analysis 
whose model dependence has not been properly identified 
and stated. On the experimental side, extracting a cen- 
troid energy necessitates the identification of all many- 
body states with a given J'" from both one-nucleon strip- 
ping and pick-up reactions, which is not often possible. 
This is particularly critical as one moves away from dou- 
bly closed shell nuclei. 

Theoretically, various levels of model dependence arise 
in the computation of ESPEs. On the deepest level, it is 
essential to understand that ESPEs depend, contrary to 
true observables, on the resolution scale A used to define 
and solve the nuclear many-body problem. As a result, 
changing A through, e.g., a unitary transformation on 
Fock space, changes ESPEs while leaving actual observ- 
ables invariant. In this sense ESPEs are similar to spec- 
troscopic factors; i.e. they can be used as a A-dependent 
analysis tool but cannot be seen as fundamental observ- 
able quantities. Moreover, and on a less fundamental 
level, approximations are often introduced in the compu- 
tation of ESPEs that generate an artificial dependence 
on the single-particle basis used. These various points 
will be discussed and illustrated in the present paper. 

Difficulties may also arise when comparing ESPEs 
computed from an ab-initio approach on the one hand 
and from more effective methods, e.g. shell model and 
energy density functional, on the other. For instance, 
while the empirical shell-model " anchors" ESPEs on one- 
nucleon addition (removal) energies to (from) the closed- 
shell core nucleus of reference, this is not the case in an 
ab-initio context, as will be illustrated below. 

The present paper follows the approach by Baranger as 
a way to delve further into the meaning and the useful- 
ness of ESPEs by addressing questions 3 and 5 above, as 
well as by quantifying the error made when using approx- 
imations to Baranger's definition. The paper is organized 
as follows. Section |lT] collects essentially known results 
regarding the definition and the computation of ESPEs. 
Such a rather exhaustive introductory part is needed to 



discuss points that have often been overlooked over the 
years. Section IIIII details the computation of ESPEs 
within the frame of the coupled-cluster (CC) method. 
Section HV] reports our results and illustrates various key 
properties of ESPEs. Specifically, the effect of correla- 
tions on both one-nucleon separation energies and on ES- 
PEs is discussed, focusing first on a few specific examples 
before addressing systematics in oxygen and calcium iso- 
topes. Starting from the ab-initio perspective provided 
by our results, the textbook rationale behind the trun- 
cated shell model is then briefly justified. Next, errors 
made by computing ESPEs in approximate ways are ad- 
dressed before illustrating the deeper model dependence 
of ESPEs associated with their intrinsic resolution scale 
dependence. Conclusions are given in Sec. El 



II. EFFECTIVE SINGLE-PARTICLE ENERGIES 

In low-energy nuclear structure theory, one usually 
starts from an independent-particle model to convey the 
basic notions of single-particle states and shell structure. 
In this context, one resorts to systems that can be pos- 
tulated a priori as being little influenced by correlations 
such that an effective independent-particle picture can be 
safely used. In a second step, actual correlations are in- 
troduced to explain, e.g., the fragmentation of the single- 
particle strength visible in one-nucleon transfer reactions. 
Such a pedagogical presentation makes it difficult to pic- 
ture the possibility to define and extract a posteriori an 
effective^ underlying single-particle shell structure in the 
presence of correlations, i.e. for A-body systems that are, 
strictly speaking, always correlated. It is thus more in- 
structive to start from a realistic picture of the nucleus, 
i.e. a rather strongly correlated system, and extract a 
posteriori an effective single-particle shell structure from 
which correlations are to a large extent, but not entirely, 
screened out [l^ . 

To do so, we introduce the nuclear Hamiltonian under 
the formi H = T + V^^ + F^n ^^ere T denotes 

the kinetic energy operator while V^^ corresponds to a 
B-body interaction. We limit ourselves to 2N and 3N 
interactions throughout the formal part of the paper and 
to 2N forces in actual applications. Studying the impact 
of 3N interactions and forces of higher rank is postponed 
to future works. Given H, eigenstates and eigenenergies 
of the A-nucleon system are obtained by solving 

H\^^) ^ E^\^^) , (1) 



^ The complication associated with the self-bound character of 
the nucleus, i.e. the need to subtract the centcr-of-mass mo- 
tion in order to deal with internal many-body states and eigen- 
energies llTll . is overlooked in the present paper. Dealing with 
this difficulty in actual calculations is mandatory but would un- 
necessarily complicate the analytical expressions presented here 
without modifying significantly the outcome. 
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where the symmetry quantum number denoting the par- 
ticle number has been singled out. The label n collects 
a principal quantum number n^, total angular momen- 
tum J^, the projection of the latter along the z axis M^, 
parity and isospin projection along the z axis of 
the many-body state of interest. Use of the Greek label 
will be made to denote the subset of quantum num- 
bers = (n^, J^,T^). Due to rotational invariance of 
the nuclear Hamiltonian, eigenenergies = ^ are 
independent of M^. 

In the following, we consider a spherical single- 
particle basis {a|,} appropriate to discussing the spher- 
ical shell structure. Basis states are labelled hy p = 
{upjiTp, jp,mp,Tp} = {np^nip^ap}, where Up represents 
the principal quantum number, iTp the parity, jp the total 
angular momentum, nip its projection along the z-axis, 
and Tp the isospin projection along the same axis. 

We also consider the direct-product basis {bt^^}, 
where f is the position vector, a the projection of the 
nucleon spin along the z axis, and r its isospin projec- 
tion. 



of motion given by jlSf 

[h°^ + E(c.)]^ 



(4) 



and similarly for {Vi,,Ej^), where (observable) one- 
nucleon separation energies arc defined through 



^ E^' 



1 

- Eo 



E' 



= Eq - E^ 



(5a) 
(5b) 



The energy- depeniieni potential S(w) denotes the dynam- 
ical part of the irreducible self-energy [31 that naturally 
arises in self-consistent Green's-function theory and that 
is to be evaluated at the eigensolution in Eq. (|4|). 
The static field is defined in Eq. ([T8| and contains 
both the kinetic energy and the cncvgy- independent part 
of the one nucleon self-energy. One can show from Eq. Q 
that the long-distance behaviour of the radial part of the 
overlap function is governed by the corresponding one- 
nucleon separation energy, e.g. for iJ+ < 



A. Spectroscopic amplitudes 

The physical processes providing information on the 
single-particle shell structure are one-nucleon transfer re- 
actions. Although the discussion can be carried out for 
the transfer on any initial [l3| . many-body state, we 
restrict ourselves in the following to the transfer on the 
ground state l^*^) of an even-even system, i.e. a J'^ = 0+ 
state. Furthermore, we consider this nucleus to be of dou- 
bly closed-shell character^. 

In this context, let us introduce {Vv) as the ampli- 
tude to reach a specific eigenstate of the 
A-l-1 (A-1) system by adding (removing) a nucleon in a 
specific single-particle state to (from) the ground state of 
the A-body system I^'q )• Such spectroscopic amplitudes 
can be defined through their representation in any given 
single-particle basis. In basis {al}, they read 



(2a) 
(2b) 



whereas their representation in basis {&|?^^} provides the 
associated wave functions or overlap Junctions 



V,y{raT) 



(3a) 
(3b) 



An important property regarding the asymptotic be- 
haviour of overlap functions derives from their equation 



Uf_,{raT) 



At 



>-\-oo 



(6) 



where denotes the so-called asymptotic normalization 
coefficient (ANC) while the decay constant is given by 
<r+ = (— 2m£'+//i^)^/^, where m is the nucleon mass^^. 
A similar result can, of course, be obtained for V^^rar) 
whose decay constant relates to E~ . 

From spectroscopic amplitudes one defines addition 5*+ 
and removal S~ spectroscopic probability matrices asso- 
ciated with states l^*^^^) and j'l'^"^), respectively. Their 
matrix elements read in basis {aj^} 



c+pq 



9 

y-P* yq 



(7a) 
(7b) 



such that their diagonal parts, when expressed in the co- 
ordinate space basis, are nothing but transition densities 
for the one-nucleon transfer from 

l*o) to 11';^+^) and 

I'l'^"-^), respectively. 

Tracing the two spectroscopic probability matrices 
over the one-body Hilbert space Hi gives access to spec- 
troscopic factors 



SF+ ^ J2 \K\' = J2 [dr\U,{rarf , (8a) 
SF- EE J2 iKfl'^E JdrmraT)f , (8b) 



^ Such a notion relates to the filling of shells in the uncorrelated, 
e.g. Hartree-Fock, picture. 



^ Subtracting the center-of-mass motion would lead to using the 
reduced mass of the added/removed nucleon. 
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which are nothing but the (basis-independent) norm of 
spectroscopic amphtudes. A spectroscopic factor char- 
acterizes to what extent an eigenstate of the A-l-1 (A-1) 
system can be described as a nucleon added to (removed 
from) a single-particle state on top of the ground-state 
of the A-nuclcon system. Such a feature intrinsically de- 
pends on the resolution scale A characterizing the nuclear 
Hamiltonian and is thus not, strictly speaking, observ- 
able [I9i, 2M- Still, spectroscopic factors can serve as a 
tool to analyse the results obtained at a given resolution 
scale. 



B. Spectral function and strength distribution 

We now gather the complete spectroscopic informa- 
tion associated with onc-nucleon addition and removal 
processes in the so-called spectral function S(a;) = 
S~'"(a;)-|-S~(w). The spectral function denotes an energy- 
dependent matrix over Tii whose elements in basis {a^} 
are defined through 

where the first (second) sum is restricted to eigenstates of 
H in the Hilbert space Ha+i CHa-i) associated with the 
A+1 (A-1) system. It is of interest to introduce the n*"^ 
moment of the spectral function that defines an encrgy- 
independent matrix over Hi through 



-\-QO 



Oj''S{uj)dLO. 



(9) 



One can easily obtain from {ap,aj} = 6pq that the 
zeroth-moment is nothing but the identity matrix 



E 



q+pq 



(10) 



such that the diagonal matrix element of S(a;) possesses 
the meaning of a probability distribution function (PDF) 
in the statistical sense, i.e. the combined probability of 
adding and removing a nucleon to/from a specific single- 
particle basis state \p) integrates to 1 when summing over 
all final states in the A±l systems. 

Last, but not least, we introduce the spectral strength 
distribution (SDD) as the trace of the spectral function 
matrix 

Siio) EE Tvn, [S(c.)] (11) 
= ^ SF+ 6{uj -E+)+J2 SF- S{u: - E~) , 

fiGHA+i v^Ha-i 

which is a basis-independent function of the energy. 

C. Independent-particle vs correlated systems 

It is of pedagogical interest to discuss the typical pat- 
terns displayed by the spectral strength distribution for 
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FIG. 1: (Color online) Schematic representation of one- 
nucleon addition and removal spectroscopic information for 
an independent-particle system. Left: binding energy for 
the ground-state of an even-even system and for the states 
of neighbouring nuclei reached by direct one-nucleon addi- 
tion and removal processes. Right: corresponding spectral 
strength distribution. 



both independent-particle and correlated systems. The 
goal of this exercise is to illustrate in what sense ob- 
servable one-nucleon separation energies cannot be inter- 
preted as single-particle energies as soon as correlations 
are present in the system. 

Figure [1] provides a schematic display of one-nucleon 
addition and removal spectroscopic information for an 
independent-particle system. As many-body eigenstates 
of H take the form of Slater determinants in such a 
case, there exists a particular single-particle basis of Hi 
in which addition and removal spectroscopic probability 
matrices read 



Q+pq — A A A 
^ fj, PQ pf^ J 

c-pq — AAA. 

^1/ ^pi/ ^p2 5 



(12a) 



where i and a characterize occupied ("hole") and unoc- 
cupied ("particle") states in the Slater determinant as- 
sociated with the A-nucleon ground-state, respectively. 
Consequently, the many-body states reached by direct 
one-nucleon addition and removal processes are in one-to- 
one correspondence with single-particle basis states. As 
a result of such a bijection, one-nucleon separation ener- 
gies are good candidates to play the role of single-particle 
energies. As a matter of fact, one has E"^ = 6afi and 
i?^ = CiSiv, where Ca and ti denote eigenvalues of the 
one-body Hamiltonian governing the uncorrelated system 
associated with unoccupied and occupied single-particle 
states, respectively. Because the SDD integrates to the 
dimension of Hi by construction, spectroscopic factors 
of the corresponding states are equal to 1 whereas they 
are equal to zero for all the remaining states that are not 
reached by the direct one-nucleon transfer. 

Let us now move to a correlated system. In any single- 
particle basis {at} of Hi, S+p"^ {S-P"^) is now different 
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FIG. 2: (Color online) Same as Fig.[T]for a correlated system. 



from zero for any combination'* of fj.,p and q {v-iP and 
q) indices. The SDD is thus fragmented as schemat- 
ically displayed in Figure [21 i.e. a larger number of 
many-body states are reached through the direct addition 
and removal of a nuclcon compared to the uncorrelatcd 
case^. Consequently, the number of peaks with non-zero 
strength in the SDD is greater than the dimension of Hi, 
which forbids the establishment of a bijection between 
this set of peaks and any basis of Hi. Accordingly, and 
because the SDD still integrates to the dimension of T-Li 
by construction (see Eq. (fTO])). spectroscopic factors are 
smaller than one. The impossibility to realize such a bi- 
jection constitutes the most direct and intuitive way to 
understand why observable one-nucleon separation ener- 
gies cannot be rigorously associated with single-particle 
energies when correlations are present in the system, i.e. 
as soon as many-body eigenstates of H differ from Slater 
determinants. 



D. Effective single- particle energies 

The discussion provided above underlines the fact that 
a rigorous definition of ESPEs is yet to be provided in 
the realistic context of correlated many-nucleon systems. 
A key question is: how can one extract a set of single- 
particle energy levels that (i) are in one-to-one correspon- 
dence with a basis of T-Li, (ii) are independent of the par- 
ticular single-particle basis one is working with, (iii) are 
computable only using quantities coming out of the corre- 
lated A-body Schrodinger equation and that (iv) reduce 
to HF single-particle energies in the HE approximation 
to the A-body problem. 

Let us make the hypothesis that ideal one-nucleon pick- 



up and stripping reactions have been performed such that 
separation energies (E'^^E~) and spectroscopic ampli- 
tudes (overlap functions) {U^(faT),Vv{raT)) have been 
extracted consistently, i.e. in a way that is consistent 
with the chosen nuclear Hamiltonian H{A) defined at a 
resolution scale A. In such a context, a meaningful defi- 
nition of ESPEs does exist and goes back to French [ll| 
and Baranger [l2|. It involves the computation of the 
so-called centroid matrix which, in an arbitrary spherical 
basis of "Hi reads 



S+P1E+ 



(13a) 



and is nothing but the first moment M^^^ of the spectral 
ftmction matrix (see Eq. [3]). Effective single-particle en- 
ergies and associated states are extracted, respectively, 
as eigenvalues and eigenvectors of i.e. by solving 



L cent ^/,ccnt 



„ccnt ^/,ccnt 



(14) 



where the resulting spherical basis is denoted as {cj^}. 
Written in that basis, centroid energies invoke diagonal 
spectroscopic probabilities^ 



E 



S+PPE+ 



E ^^"""^^ 



(15) 



and acquire the meaning of an average of one-nucleon sep- 
aration energies weighted by the probability to reach the 
corresponding A-l-1 (A-1) eigenstates by adding (remov- 
ing) a nucleon to (from) the single-particle state ipp°'^^. 
Centroid energies are by construction in one-to-one cor- 
respondence with states of a single-particle basis of Hi 
which, as already pointed out before, is not the case of 
correlated one-nucleon separation energies with non-zero 
spectroscopic strength. 

Equation (fTi)) ensures that ipp"'^*^{faT) and e™"' are 
consistent in the sense that the asymptotic behaviour of 
the former is driven by the latter, e.g. for ep°"' < the 
radial part of the wave function behaves asymptotically 
as 



^ r-> + oo r 



(16) 



where = (-2me™ /?i ) ' . Such a result under- 
lines that single-particle wave-functions associated with 
ESPEs are intrinsically different from overlap functions 
Ufi{raT) (Vi/(r(Tr)) which are probed in transfer experi- 
ments. 



* Except for selection rules dictated by symmetries tliat lead, ac- 
cording to Eq. I I20I I. to TTp = TTp, jp = and Tp = — Tq. 

^ Of course, the dimension of Ha+i or Ha—1 remains the same 
whether the system is correlated or not. 



The definition of 6^°"' sometimes incorporates the denominator 



E 



in Eq. I I15I I to compensate for 



the possibility that, e.g. experimentally, normahzation condi- 
tion [TO] might not be exhausted. 
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Experimentally, the extraction of ESPEs requires to 
collect the full spectroscopic strength up to high enough 
missing energies, i.e. the complete set of separation en- 
ergies and cross sections from both onc-nucleon stripping 
and pickup reactions. This unfortunately limits the pos- 
sibility to perform sound comparisons on a systematic 
basis. Indeed, there are at best only a few nuclei along a 
given isotopic or isotonic chain that are characterized by 
complete enough spectroscopic data. 



E. Sum rule 

It is tedious but straightforward to prove that the n"^ 
moment of S(a;) fulfils the identity 

n commutators 



M(") - «|{[...[[a„if],if],...],at}|0. (17) 



^3N 



to- 



Using the second quantized form of T, V , and V' 
gether with identities provided in Appendix |X] and sym- 
metries of interaction matrix elements, Eq. (jl7p applied 
to n = 1 leads to [il,[il,[2ll 

^cent ^ J. , f72N [1] , i y3N [2] 

"'pq -^PQ ' / J ^ prqs rsr ' ^ / ^ ^ prtqsv rsvrt 

rs rstv 

= h°° , (18) 

where Vp^^ and Vp^gsv are anti-symmetrized matrix el- 
ements and where 

^ {^^\alap\^^}^Y.^^*V^ , (19a) 

(19b) 



denote one- and two-body density matrices of the corre- 
lated A-body ground-state, respectively. The static field 
h°°, already introduced in Sec. Ill A| contains both the 
kinetic energy and the energy-independent part of the 



onc-nucleon self-energy in the A-body ground state [21 
Equation ()18|) demonstrates that the centroid matrix 



is a one-body field possessing a simple structure and an 
intuitive meaning. In particular, the centroid field re- 
duces to the Hartree-Fock (HE) mean field in the HE 
approximation. As a result, ESPEs are nothing but HE 
single-particle energies in such a case and are equal to 
onc-nucleon separation energies according to Koopmans' 
theorem (2^. Consistently, overlap, centroid, and HE 
single- particle wave- functions coincide in that limit. Of 
course, centroid energies also reduce to eigenvalues of 
the one-body Hamiltonian in the limit of an uncorrelated 
system. When correlations beyond HE are switched on, 
ESPEs arc modified through the presence of correlated 
density matrices in Eq. (jlSp : i.e. the B-nucleon interac- 
tion is folded with the correlated (B-l)-body density ma- 
trix Through that transition, ESPEs continuously 
evolve as centroid energies rather than as observable sep- 
aration energies such that Koopmans' theorem does not 



hold any more. Centroid energies are schematically com- 
pared to observable binding and separation energies in 
Eigurc 131 

On the practical side, Eq. underlines that the av- 
eraged information contained in ESPEs only requires the 
computation of the A-body ground-state. As long as 
one is not interested in the full spectroscopic strength 
of the A±l systems but only in their centroids, one only 
needs to compute one nucleus instead of three. In prac- 
tice however, Eq. ((22)) is rarely computed in terms of 
the correlated density matrix, e.g. shell-model applica- 
tions usually invoke a filling approximation typical of an 
independent-particle approximation. This is believed to 
be a decent approximation as long as (i) low-lying states 
carry a major part of the single-particle spectroscopic 
strength, as for the transfer on a doubly closed-shell nu- 
cleus, and (ii) nucleons of the other species are themselves 
not strongly correlated, because of pairing for example. 
See, e.g., Ref. [2^ and references therein for a related 
discussion. Such an issue becomes critical whenever one 
is looking into, e.g., the neutron shell structure of a neu- 
tron open-shell nucleus. In such a situation, a normal 
filling is inappropriate and it is mandatory to fold the 
monopole interaction in Eq. (j22p with a density matrix 
reflecting the presence of correlations in the system. 

Using that the even-even ground state the one-nucleon 
transfer is performed on is a j'^ = 0+ state, Wigner- 
Eckart's theorem allows one to obtain the explicit de- 
pendence of spectroscopic amplitudes on rup and M^, 
i.e. 

UP EE t/„";["''U.^„^<5M.™,, (20a) 
yp ^ vZ-^'^-U^.^ajAU-mA-^r^^ (20b) 

such that the single-particle operator picks out the angu- 
lar momentum, the parity and the isospin projection of 
the A±l state the transfer goes to; i.e. jp = J^, TTp = H^ 
and Tp = Tfj^ — Tq. Consequently, the one-body density 
matrix of the A-body ground state reads 

Fpq — yiipUq "apCfg ^Jmprrig , 



(21) 



such that, retaining the 2N force only for simplicity and 
expressing its anti-symmetrized matrix elements in a jj- 
coupled scheme, one obtains [l2| - [l3 | in basis {c|,} 



flap 



EV^ -[apQgQpa^] [Qg] 



(22) 



where is the reduction of the 2N interaction 

to its so-called monopole, i.e. angular averaged, part. 
Higher multipoles and in particular the quadrupole part 
that drives the dominant part of correlations are screened 
out from ESPEs. 



F. Resolution-scale dependence 

Let us briefly explain the intrinsic resolution-scale de- 
pendence of ESPEs. Such a feature derives directly 
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FIG. 3; (Color online) Schematic picture. Top: total bind- 
ing energies (Eq. [l]) of three successive nuclei and associ- 
ated one-nucleon addition/removal energies (Eq. [5} from the 
ground state of the intermediate system. Bottom left: spec- 
tral strength distribution (Eq.[TT]). Bottom right: correspond- 
ing ESPE spectrum (Eq. [15}. The color coding underlines 
that ESPEs close to the Fermi energy contain significant con- 
tributions from both addition and removal channels. 



from the resolution scale of spectroscopic amplitudes |l9l . 
[20| entering the definition of the centroid matrix, see 
Eq. ()13ap . Following the philosophy of the similarity 
renormalization group (SRG) [2J, [23, we consider a 
change of resolution scale via a unitary transformation 
U (A) of the Hamiltonian 



H{A) = U{A)HU\A) 



T + V^^ {A) + V'' (A) + 



3N( 



(23a) 
(23b) 



such that^ the associated observable, i.e. the eigenen- 
ergy, remains unchanged E-^{A) = E^. Similarly, any 
observable associated with a Hermitian operator O must 
remain invariant, which imposes the transformation of O 
according to 0(A) = U{A) O C/t(A). 

Let us now come to ESPEs. The key difference from 
an observable resides in the fact that the very nature of 
ESPEs is to inform us of effective single-nucleon degrees 
of freedom inside the nuclear medium, independently of 
the form of the Hamiltonian. In other words, the choice 
is made to keep the definition of ESPEs independent of 
A. Before or after transformation [23al ESPEs arc always 
extracted through Eq. (|13ap . where 5'^p'?(A) retains the 
same formal expression as before, i.e. they invoke spec- 
troscopic amplitudes computed through 



K?(A) 



(v&^+i(A)|at|*^(A))* 
(vI'A-i(A)|a,K(A))*. 



(25a) 
(25b) 



Contrary to the many-body states involved, operators aj, 
and Op are not transformed in Eq. p5|) . which generates 
automatically an intrinsic dependence of and on 
A. One could, of course, choose to transform operators 
and Up in the definition of spectroscopic amplitudes in or- 
der to make the latter invariant under the unitary trans- 
formation. However, transforming a^, e.g., would result 
into a linear combination of operators of the form aj, 
aJaJ^Os, . . . such that the spectroscopic amplitude would 
not provide the information one was after in the first 
place, e.g. the overlap between eigenstates of H{A) in 
the A-l-1 (A-1) system and the state obtained by adding 
(removing) a nucleon to (from) a given single-particle 
state \p) on top of the A-body ground-state. Eventually, 
the resolution-scale dependence of spectroscopic ampli- 
tudcs propagates to their norm, i.e. spectroscopic fac- 
tors [13, [2I, and to ESPEs. 

The discussion provided above points to an important 
conclusion. The information one is sometimes after, e.g. 
computing spectroscopic factors and ESPEs, is not nec- 
essarily observable. Such an information is not abso- 
lute and can be modified by a redefinition of inaccessible 
quantities, i.e. the Hamiltonian and its eigenvectors in 
the present case, which leaves of course true observables 
untouched. It remains to be seen how much ESPEs are 
changed in actual calculations by varying the resolution 
scale A over a reasonable interval of interest. This ques- 
tion is addressed in Sec. lIVI It could very well be that the 
induced variation of the ESPEs is negligible compared 



where the scale characterizing the initial H is omitted for 
simplicity. As can be trivially shown, Eq. ([23a|) induces 
a transformation of eigenvectors of Eq. ([1]) 



|*^(A)) ^ [/(A) \^^) , (24) 



^ In practical applications, such an invariancc is broken to some 
extent due to the approximate way of performing the transfor- 
mation of the Hamiltonian, e.g. neglecting induced many-body 
interactions in Eq. I)23b|l , and to approximations performed when 
solving the A-body problem I26l l27ll . However, the discussion of 
the present section is concerned with tracking what happens in 
the hypothesis of an exact unitary transformation and an exact 
solution of the A-body Schroedinger equation. 
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to other sources of uncertainties, e.g. approximations in 
their computation. Still, it is of prime importance to keep 
such an intrinsic model dependence of ESPEs in mind. 



III. COUPLED CLUSTER METHOD 

One-nucleon separation energies and spectroscopic am- 
plitudes introduced respectively in Eq. ([5]) and Eqs. ([2][3]) 
are defined without any reference to a particular method 
used to solve the many-body problem. 

We are presently interested in using the ab-initio 
coupled-cluster method (CCM). Let us briefly outline 
the procedure to compute ground and excited states of a 
closed (sub-)shell nucleus A and of odd Ail neighbours 
within CCM. From there, all needed quantities to com- 
pute ESPEs can be extracted. In CCM, the exact ground 
state is written in the exponential form 

l*o) = e^l<i>o), (26) 

where |$o) is an uncorrelated single-reference Slater 
determinant built from a convenient spherical single- 
particle basis, usually chosen as mean-field HF orbitals. 
Many-body correlations beyond the mean field are intro- 
duced by the operator T = Ti + T2 + . . . + Ta, which is a 
linear expansion in n-particle-n-hole excitation operators 
r„, with n = 1, . . . , A. 

The only approximation occurring in CCM regards the 
truncation of T to a given low-lying excitation level; 
e.g. T w Ti + r2 is the most commonly used ap- 
proximation known as the coupled-cluster method with 
single and double excitations (CCSD). Inserting the 
coupled-cluster ansatz ([25)1 into the A-body Schrodinger 
equation (Eq. ([T])) and projecting from the left with 
{^o\e~'^ , {^f\e~'^ and {^f!j'\e~'^ respectively, coupled- 
cluster equations are obtained under the form 

($o|e-^i/e^|$o) = E^, (27a) 
($?|e-^7Ie^|$o) = 0, (27b) 
($^/|e-^iJe^|$o) - 0. (27c) 

These equations determine the unknown amplitudes en- 
tering Ti and T2 as well as the ground state energy E^. 
Here ($"| and {^fj\ are one-particle-one- hole and two- 
particle-two-hole excited reference states. 

Equation (P7|) underlines that the similarity- 
transformed Hamiltonian H = e~^He^ plays a key role 
such that its ground state is nothing but the reference 
state |$o)- The operator H is not Hermitian, which 
implies that coupled-cluster theory is manifestly non- 
variational. The non-variational nature of CCM makes 
it necessary to access both right and left eigenstates of 
H to compute associated one- and two-body density ma- 
trices. Such eigenstates of H can be computed through 
the so-called equation-of-motion coupled-cluster method 
(EOM-CCM)^ The idea of EOM-CC is essentially to 
diagonalize H within a subspace of n-particle-m-hole 



excited reference functions. Within the EOM-CCSD 
approximation, right and left eigenstates of closed-shell 
nucleus A are given by 

K) - R^l'^o) , (28a) 
{Lt\ - (*o|i;^ , (28b) 

where i?^ (L^) is a linear combination of one-particle- 
one-hole and two-particle-two-hole (de-)excitation oper- 
ators. Similarly, EOM-CC is the method of choice to ac- 
cess eigenstates of odd Ail neighbouring nuclei accord- 
ing to = R^^'\<^>^) and {L^^'\ = {^^\L^^\ 

where now R^^^ i^^^'^) denotes a linear combination 
of one-particle (one-hole) and two-particle-one-hole (one- 
particle-two-hole) (de-)excitation operators (see for ex- 
ample Ref. [2^ for further details). Left and right eigen- 
states form a bi-orthogonal set, i.e. 

(i^|i?^') = V- (29) 

where we have dropped the superscript referring to nu- 
cleus A and Ail. Right eigenstates i?^ arc solutions of 
the eigenvalue problem 

iHR^)c\^o) = E^R^l^o) , (30) 

and similarly for left eigenstates L^. Here, {HR^)c de- 
notes all terms that connect H with R^. The one- and 
two-body density matrices of the A-body ground state 
together with onc-nucleon spectroscopic amplitudes and 
probabilities can now be computed according to 



pW = ($o|L^4ap|<f>o) , (31a) 
P^pI^^ = {(^o\Laalalaqap\^o) , 

and to 

Sp'^ ^ ($o|i^a^i?^^|<f>o)($o|^^Ml*o) ,(32a) 

^ ($o|iMi?^'|$o)($o|i^'a^|*o) .(32b) 

Using the Baker-Campbell-Hausdorff commutator ex- 
pansion, one can derive finite and closed-form alge- 
braic expressions for similarity-transformed operators a^, 

and alalaqap. See Refs. [2^ [s^l for details 
on derivation and computation of spectroscopic factors 
within coupled-cluster theory. 

IV. RESULTS 

Results shown below have been obtained using a 
2N force only, thereby omitting forces of higher rank. 
In order to improve convergence properties and make 
the nuclear many-body problem more perturbative, we 
use a soft Mow-k 2N interaction [3l| obtained with a 
smooth regulator [s^ for various cutoff values between 
A =2.0 fm^^ and 3.0 fm^^. These soft Vjow-k interactions 
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are obtained by evolving down the N'^LO Chiral interac- 
tion [3^ with cutoff =500 MeV. For the single-particle 
model space, we use the HF basis built from N-|-l= 13 
major oscillator shells, with a fixed oscillator frequency 
of hid =16McV. This model space is sufficient to obtain 
fully converged results for medium mass nuclei with soft 
Mow-k interactions (see Ref.flTj). 



A. Turning on correlations in a controlled way 

Let us first illustrate in a pedagogical manner the effect 
of correlations on both one-nucleon separation energies 
and ESPEs. To do so, we apply Wick's theorem with 
respect to the HF vacuum |$[^^) and write the Hamilto- 
nian in normal-ordered form using the HF single-particle 
basis {dj^}, i.e. 



H = h 



HF 



where 



h 

together with 



HF _ 



HF 



E^^ + Y^e^^ -.did,: 

p 

1 V • dUU d • 



(33) 

(34a) 
(34b) 



2N 

2 ' P1P1 

p,q=l 



,HF 



Tpp + Y.% 



2N 
pqpq 



(34c) 
(34d) 



9=1 



Scaling the residual interaction V^cs by a factor A £ [0, 1], 
one defines a parameter-dependent Hamiltonian Hx = 
/i^^+A Vtcs that tunes correlations between the two limits 
of interest, i.e. from the uncorrelated regime Hq = 
to the fully correlated regime Hi = H . Eventually, we 
solve EOM-CC equations repeatedly for several values of 
A S [0, 1] and for the specific cutoff value A =2.4 fm~^. 

As a first example. Figure S] displays, as a function 
of the residual interaction strength, one-neutron sepa- 
ration energies between ^^O ground-state and low-lying 
states in ^^O along with corresponding ESPEs. Plot- 
ted separation energies correspond to the (main) lowest 
peak in the additional sector of the SDD for each 
symmetry block. As expected, one-neutron separation 
energies and ESPEs are equal in the uncorrelated limit 
(A = 0) and are nothing but HF single-particle energies, 
i.e. Koopmans' theorem is fulfilled. Turning on correla- 
tions, two important features manifest themselves. First, 
the SDD is fragmented such that the separation energy 
of the state carrying the largest strength for a given 
goes down significantly. Second, correlations only slightly 
impact centroid energies that keep a strong memory of 
HF single-particle energies. Thus, although ESPEs are 
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FIG. 4: (Color online) One-neutron addition energies on 
^^0 ground-state and corresponding ESPEs Bp''"* as a function 
of the residual interaction strength. 

1 
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FIG. 5: (Color online) Removal spectroscopic factor SF , 

1/2 

of the lowest 1/2^ state in ^'^O as a function of the residual 
interaction strength. 



not independent of correlations, the latter are essentially 
screened out as discussed earlier. Eventually, one-nucleon 
separation energies and corresponding ESPEs can differ 
by several MeVs. This clearly points to the fact that sep- 
aration energies should not be identified as ESPEs and 
vice versa. 

Let us now include more detail by focusing on the low- 
est 1/2+ state in ^^O, which can be accessed by remov- 
ing a neutron from the ground state of '^'^O. Figure [5] 
displays the corresponding spectroscopic factor SF~^^_^ 
as a function of the residual interaction strength. While 
^^1/2+ = 1 for A = as expected, it decreases gently as 
the residual interaction is switched on to reach a value 
of 0.92 in the fully correlated case. Using it as a (scale- 
dependent) analysis tool, such a spectroscopic factor tells 
us that the lowest 1/2+ state keeps a well-pronounced 
single-particle character even in the fully correlated limit. 

Figure [S] shows corresponding neutron separation en- 
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FIG. 6: (Color online) Same as Figure [5] for one-neutron re- 
moval energy E~i^^ and ESPE ea^'a- 




22 24 
Oxygen Isotope 

FIG. 7: (Color online) Evolution of selected one-neutron sep- 
aration energies Et, and corresponding ESPEs gp''"* from 
to 2«0. 



ergy and ESPE as a function of A. In the uncorrelated 
limit, i.e. A = 0, the separation energy and ESPE co- 
incide, while for larger A they start to deviate. While 
e2™*2+ is O'^ly slightly influenced by correlations, 
dives significantly as A increases from to 1. Eventually, 
correlations add about 1.5 MeV to the separation energy, 
such that it differs from e^7/*2+ 1.7 McV for A = 1. 
Even though the 1/2+ state retains to a large extent 
its single-particle nature, its energy is strongly impacted 
by correlations and does not provide clean information 
about the effective single-particle shell structure. 

B. Systematics in oxygen and calcium isotopes 

We now discuss the evolution and trends of low-lying 
one-neutron addition and removal energies together with 
ESPEs in doubly closed-shell oxygen and calcium iso- 
topes. The present calculations are performed with the 
specific cutoff value A = 2.4 fm~^. As a result, the neu- 
tron drip line in oxygen and calcium isotopes is wrongly 
predicted to be located beyond and ^°Ca, respec- 
tively. Three-body forces seem to be mandatory to cor- 
rectly reproduce the drip line location at ^^O [1] for oxy- 
gen isotopes. It remains to be seen in which way forces 
of higher rank modify qualitatively or quantitatively the 
conclusions of the present investigation. 

Figures [7] and [5] show that separation energies 
(£"+, E~) systematically and significantly differ from cor- 
responding centroid energies. As for the energetics, these 
results illustrate that ab-initio approaches describe dou- 
bly magic nuclei such as ^^'^^o and ^^'''^Ca as strongly 
correlated systems. Most importantly, how much sepa- 
ration energies differ from centroid energies significantly 
depends on the nucleus/state, in a way that cannot eas- 
ily be traced back to one specific feature. Consequently, 
opening or closing of shell gaps in the separation energy 
spectrum are not in one-to-one relationship with those 
emerging in the ESPE spectrum. One does observe that 




Calcium Isotope (A) 



FIG. 8: (Color online) Same as Figure [7] from ''°Ca to "^"Ca. 



^^O and °Ca display the strongest correlations of all, 
which may be related to their N=Z character. Trac- 
ing the isospin dependence of correlations in oxygen and 
calcium isotopic chains, there seems to be a systematic 
trend with increasing asymmetry N-Z. In both oxygen 
and calcium chains, correlations become less important 
for the neutrons close to the Fermi surface when increas- 
ing isospin asymmetry N-Z. This trend is consistent with 
Ref. [3J|, where it was found that the spectroscopic fac- 
tor for removing a neutron close to the Fermi surface in- 
creases with increasing isospin asymmetry and is close to 
one for ^®0, while the spectroscopic factor for removing 
the outermost protons is largely quenched with increas- 
ing isospin asymmetry. 

One can conclude from Figs. [7] and [5] that inferring one- 
nucleon separation energies from ESPEs is not straight- 
forward, even in doubly closed-shell nuclei. One should 
thus simply not use one for the other. 
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C. Effective shell model 

Equation 1181 was obtained following an ab-initio strat- 
egy, i.e. considering all nucleons as active and interact- 
ing via realistic 2N and 3N interactions in a large enough 
single-particle Hilbert space. In the traditional effective 
shell-model, however, the equivalent of Eq. is derived 
from an effective Hamiltonian defined for n™' active nu- 
cleons in a restricted valence space above a closed core 
composed of Tf-™'' nucleons and below an excluded space. 

In such a context, ESPEs are "anchored" on (exper- 
imental) one-nucleon addition energies to the core nu- 
cleus, i.e. 6™"^° = i?+(5pfc. As seen in previous sections, 
this constitutes in fact a bad approximation when tak- 
ing an ab-initio perspective. However, and as confirmed 
by the present investigation, low-lying states obtained by 
adding (removing) one nucleon to (from) a doubly closed 
shell nucleus do possess a well-defined single-particle 
character. The reason why the corresponding ESPE dif- 
fers significantly from the separation energy is due to the 
fact that the former collects small strength rejected to 
rather high missing energies. Eventually, the fact that 
low-lying states carry most of the strength makes them 
good candidates to represent quasi single-particle degrees 
of freedom. Such a textbook result constitutes the ba- 
sic justification for the effective shell model that omits 
the high-lying fragmented strength and recollects the full 
strength into low-lying states of the core-|-one-nucleon 
system. 



D. Using a fixed single- particle basis 

Computing ESPEs in an approximate fashion gener- 
ates a model dependence that may compromise their use- 
fulness. It is, for example, customary to use uncorrelated 
occupations of single-particle states in place of the cor- 
related one-body density matrix in Eq. ((T5)) and/or to 
define ESPEs as the diagonal matrix elements of the cen- 
troid field /i"^""* in an a priori chosen single-particle basis, 
e.g. a harmonic oscillator basis, rather than as its eigen- 
values. The latter approximation is formally questionable 
as it provides a basis-dependent definition of ESPEs, the 
quality of which depends on the realistic character of the 
chosen basis. In practice, the quantitative impact of such 
an approximate scheme depends on the situation. 

Figure IHl compares in oxygen isotopes properly com- 
puted ESPEs with diagonal matrix elements of in 
the HF basis used in the calculation® . As can be inferred 



A more drastic approximation not shown here consists in using 
diagonal matrix elements in a harmonic oscillator basis. This is 
the choice usually made within the frame of the interacting shell 
model. In practice the model space is usually small enough to 
contain only one state per symmetry block, e.g. {l,j,m) block in 
spherical symmetry. In such a case the arbitrarily chosen basis 
is necessarily the eigenbasis of the ccntroid field, underpinning 



from the comparison with Fig. [71 the approximation in- 
duces errors on ESPEs that are of the same order as their 
difference with one-nucleon separation energies and that 
arc, in some cases, significant relative to their absolute 
values. The error depends both on the state and on the 
system, i.e. it might go in opposite directions depending 
on the state and/or the nucleus under consideration. 

Interestingly, there exists cases for which the ordering 
of approximate ESPEs at the Fermi level is inverted com- 
pared to full-fledged ones, e.g. for 2si/2 and 1^5/2 levels 
in ^^O. Knowing that full-blown ESPEs reproduce the 
ordering of one-neutron separation energies across the 
whole set of oxygen and calcium isotopes, such an inver- 
sion is of noticeable importance. The inversion seen in 
2^0 is consistent with Fig. 2a of Ref. H where ESPEs 
were computed as diagonal matrix elements of /i'^*'"* in an 
a priori chosen harmonic oscillator basis. The fact that 
2si/2 and 1^5/2 levels are actually not inverted in Fig. 
2a of Ref. Q is simply due to the fact that ESPEs are 
anchored on empirical values in ^^O. Correcting mentally 
for such a fact, one recovers the level inversion seen for 
"Diag-ESPE" in Fig. Oof the present paper. 




Oxygen Isotope 

FIG. 9: (Color online) ESPEs e^°°' compared to diagonal 
matrix elements of the centroid field ft,™' in the underlying 
HF basis. Results are displayed from ^''O to ^^^O. 



E. Resolution scale dependence 

A more fundamental model dependence of ESPEs that 
remains even when computing them as eigenvalues of the 
centroid matrix relates to the resolution scale character- 
izing the Hamiltonian. We start from a Chiral Hamilto- 
nian built with a cutoff (e.g. 500 MeV here) up to 
a given order (e.g. N'^LO here). This in itself carries a 
truncation error with respect to using the complete EFT 



the strong impact of using a severely restricted model space. 
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Lagrangian. Still, this constitutes our reference Hamil- 
tonian, which at N'^LO contains both 2N and 3N inter- 
actions. In a second step, the resolution scale of the 
Hamiltonian is lowered to a value A through a renor- 
malization group transformation, defining in this way 
H{A) = -ffiow-k- Doing so softens the interactions and in- 
duces multi-body forces, e.g. 3N interactions are induced 
from the original 2N one, while preserving the original 
truncation error. As A is lowered, true observables re- 
main the same as with the original Chiral Hamiltonian 
as long as induced interactions are kept in the calculation 
and the many-body problem is solved exactly. Contrarily, 
even in such conditions non-observables quantities such 
as ESPEs are modified when changing A. This consti- 
tutes the intrinsic scale dependence of ESPEs discussed 
in Sec. IIIFI and that we presently wish to characterize. 
Of course, whenever induced interactions are discarded 
and/or the many- problem is not solved exactly, both ob- 
servable and non-observable quantities acquire an addi- 
tional artificial dependence on A. 

As original and induced three-body forces, as well as 
clusters beyond singles and doubles, are discarded in the 
present calculation, ESPEs display the two sources of 
A dependence. In order to extract the intrinsic one, one 
must first pin down the artificial scale dependence to sub- 
tract it eventually. By definition, the latter can be ac- 
cessed by focusing on true observables. Figure [10] displays 
one-neutron removal energies with J'^ = l/2+,5/2+ in 
^^O for various values^ of the momentum cutoff A of 
the 2N interaction Vjow-k- Lowering A from 3.0 to 2.0 
fm""'^ changes one-neutron removal energies by about 7 
MeV. Eventually, including induced many-body interac- 
tions, i.e. three- and possibly four-body forces p6l |27|. 
and including triples will remove such an artificial depen- 
dence of one-neutron removal energies on A. 

Figure [TO] also shows ESPEs e|™*, and e^X^^^ in 2*0. 
Clearly, they display a significantly larger cutoff variation 
than corresponding one-neutron removal energies. Such 
a feature, visible in all isotopes and for all states, is iden- 
tified with the additional intrinsic scale dependence of 
ESPEs. Mentally subtracting the cutoff dependence of 
one-neutron removal energies, one sees that such an in- 
trinsic scale dependence increases with A as the system 
becomes less and less perturbative, making ESPEs differ 
more and more from separation energies. Quantitatively 
speaking, the intrinsic cutoff dependence of e|™*^ and 
^1^5/2 ^''^ou'^ts to about 6 MeV when varying A from 2.0 
to 3.0 fm^^, which is obviously significant. More specifi- 



We keep the oscillator frequency fixed at huj =16 MeV in present 
calculations. For large cutoff values, the optimal oscillator fre- 
quency is larger, e.g. huj =24McV for A = 2.()fm^^, such that 
corresponding values shown in Fig. 1101 are not fully converged, 
e.g. they changed by about 100 keV for A = 2.6 fm~^ when using 
huj =24 MeV. Conclusions of the present section are however not 
modified by using fully converged values. 
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FIG. 10; (Color online) Neutron ESPEs and removal energies 
in ^^O for J'" = l/2+,5/2+ and A £ [2.0,3.0] frn'^ 

cally, one notes that both ESPEs do not vary identically 
across the range of A values. As a matter of fact, one 
observes an inversion of the ESPE ordering that is not 
reflected in one-neutron separation energies. Besides re- 
moving (most of) their artificial cutoff dependence, it will 
be of interest to see how much 3N forces modify the in- 
trinsic scale dependence of ESPEs. It is anyway likely 
that the latter will remain significant. 

The above result demonstrates that ESPEs are not 
absolute and can be changed significantly by modifying 
mildly the character of the Hamiltonian, i.e. by varying 
A over a rather limited range of values, while keeping 
true observables invariant. Consequently, extracting the 
single-particle shell structure and its evolution, e.g. with 
isospin, from experimental data is an illusory objective. 
However, it remains possible to perform a meaningful, 
i.e. internally consistent, analysis of a set of experimental 
data by extracting ESPEs through consistent structure 
and reaction models based on the same nuclear Hamil- 
tonian. Eventually, conclusions regarding the extracted 
shell structure, e.g. its evolution with isospin, will how- 
ever necessarily remain resolution-scale dependent. In 
particular, working at large scales is somewhat inappro- 
priate in the sense that the extracted shell structure will 
not fit with the phenomenological low-energy picture in 
such a case [2l[. 

V. CONCLUSIONS 

The present work discusses, from an ab-initio stand- 
point, the definition, the meaning, and the useful- 
ness of ESPEs in doubly closed shell medium-mass nu- 
clei. Illustrating the various points with state-of-the- 
art coupled-cluster calculations, the following conclusions 
are reached. 

• A meaningful single-particle shell structure fulfill- 
ing a minimal set of properties and known lim- 
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its, such as being independent of the particular 
single-particle basis one is working with, can be 
extracted from correlated one-nucleon separation 
energies and associated spectroscopic amplitudes. 
Such a definition relates effective single-particle en- 
ergies (ESPEs) to the so-called centroid eigenvalues 
introduced by Baranger [T2| . 

The corresponding non-interacting problem is gov- 
erned by the one-body centroid field /i'^''"*, which 
sums the kinetic energy and the energy-independent 
part of the irreducible one-nucleon self-energy that 
naturally arises in self-consistent Green's-function 
methods. 

It is customary in low-energy nuclear theory to 
compute ESPEs in an approximate fashion, e.g. by 
defining them as diagonal matrix elements of 
in an a priori chosen single-particle basis rather 
than as its eigenvalues. We have illustrated the 
fact that such approximations are unsafe. 

Even when fulfilling the required set of minimal 
properties, ESPEs are not strictly observable as 
they intrinsically depend on the resolution scale 
A of the Hamiltonian, i.e. they change under a 
unitary transformation of the Hamiltonian while 
true observables remain invariant. We have indeed 
demonstrated that ESPEs vary substantially when 
modifying mildly the resolution scale, i.e. when 
scanning a rather limited range of A values while 
correcting for the artificial dependence due to the 
omission of induced short-range many-forces. Such 
a result demonstrates that the objective of extract- 
ing a unique single-nucleon shell structure from cor- 
related observables, e.g. pinning down the nuclear 
shell evolution from experimental data, is intrin- 
sically illusory. Still, it is possible to perform a 
consistent analysis of experimental data and ex- 
tract a meaningful shell structure. To do so, one 
must use consistent structure and reaction models 
based on the same nuclear Hamiltonian. Eventu- 
ally, conclusions regarding the extracted shell struc- 
ture will anyway remain resolution-scale depen- 
dent; i.e. two practitioners using (the same) con- 
sistent method but starting from different, though 
unitarily equivalent, Hamiltonians will extract dif- 
ferent single-nucleon shell structures from identical 
observables, e.g. spectra and cross sections. This 
constitutes a puzzling but important result that 
sheds a new light on how one should look at the 
single-particle shell structure. In particular, work- 
ing at large scales is somewhat inappropriate in the 
sense that the extracted shell structure will not fit 
with the phenomenological low-energy picture in 
such a case. 

Extracting an effective single-particle shell struc- 
ture is often done for interpretation/analysis pur- 
poses and sometimes done to infer the behaviour of 



actual observables that are believed to be strongly 
correlated to patterns in the ESPE spectrum. In 
the present paper, we have focused on one-nucleon 
separation energies to low-lying states around good 
closed-shell nuclei. The conclusion is that correla- 
tions are too strong, even with low-scale interac- 
tions, for such separation energies to be in quanti- 
tative (sometimes qualitative) correspondence with 
effective single-particle energies around the Fermi 
energy. This is true even for states that retain a 
strong single-particle character, i.e. states carrying 
spectroscopic factors close to one. In a forthcoming 
study, the same type of analysis will be performed 
in connection with the energy of the 2+ excited 
state in good closed shell nuclei, i.e. we will study 
how much such excitation energies correlate with 
the Fermi gap in the ESPE spectrum. 

• The present study was conducted on the basis of 
two-nucleon interactions only. It remains to be seen 
to which extent forces of higher rank modify our 
conclusions. At the price of computing ESPEs cor- 
rectly, i.e. as eigenvalues of the centroid matrix 
rather than as its diagonal matrix elements in an 
a priori given (harmonic oscillator) basis, the au- 
thors of Refs. II] could easily repeat the present 
analysis within the frame of the shell model and 
characterize the impact of three-nucleon forces in a 
systematic way. 

• In the present work, ESPEs were defined based 
on the hypothesis that eigenstates of the nuclear 
Hamiltonian are also eigenstates of the particle 
number operator. Ab-initio calculations of open- 
shell nuclei are currently being developed on the 
basis of many-body methods breaking particle- 
number symmetry, i.e. using methods formulated 
over Fock space rather than over the Hilbcrt space 
associated with a definite number of particles. This 
is the case of the so-called self-consistent Gorkov- 
Green's function theory [s^. Extending the defi- 
nition of ESPEs accordingly [35|, one will be able 
to address properties of ESPEs in open-shell nuclei 
and conclude on their relevance in such a context. 
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Appendix A: Useful identities 

Using Wick's theorem, one can demonstrate the fol- 
lowing identities 



{[ap,alas],al} ^ +SprSqs , 
{[ap,alalatay],al} = +6pr Sq^ alat - Spr Sqt altty 
-Sps 5qy alat + 6ps 5qt alay , 

^ps ^qu Oj^^CL^CLy^dy -\- Spg 5qy Q^d^Q^ttu ^ps ^qw oi^Q^Clydu 
~^^pt ^qu Oj^j-CL^gQ^w^v ^pt ^qv ^I^Q-s^uJ^u ~1~ ^pt ^qw O'^Cl^O'v^u 
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